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ABSTRACT 

There are several supernova remnants which contain unresolved X-ray sources 
close to their centers, presumably radio-quiet neutron stars. To prove that these 
objects are indeed neutron stars, to understand the origin of their X-ray radiation, 
and to explain why they are radio-quiet, one should know their periods and period 
derivatives. We searched for pulsations of the X-ray flux from the radio-quiet 
neutron star candidate RX J0822-4300 near the center of the Puppis A supernova 
remnant observed with the ROSAT PSPC and HRI. A standard timing analysis of 
the separate PSPC and HRI data sets does not allow one to detect the periodicity 
unequivocally. However, a thorough analysis of the two observations separated by 
4.56 yr enabled us to find a statistically significant period P ~ 75.3 ms and its 
derivative P ~ 1.49 x 1CU 13 s s _1 . The corresponding characteristic parameters of 
the neutron star, age r = P/(2P) = 8.0 kyr, magnetic field B = 3.4 x 10 12 G, and 
rotational energy loss E = 1.4 X 10 37 erg s _1 , are typical for young radio pulsars. 
Since the X-ray radiation has a thermal-like spectrum, its pulsations may be due 
to a nonuniform temperature distribution over the neutron star surface caused by 
anisotropy of the heat conduction in the strongly magnetized crust. 



Subject headings: pulsars: individual (RX J0822-4300) - stars: neutron - supernovae: 
individual (Puppis A) - X-rays: stars 



Introduction 



Observations with the Einstein, ROSAT and ASCA missions have revealed several isolated, 
radio-quiet neutron star (NS) candidates in supernova remnants (see Brazier &; Johnston 1998, 
for a recent review). These objects are characterized by a lack of observable radio emission, 
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thermal-like X-ray spectra with typical (blackbody) temperatures of ~ (1 — 5) x 10 6 K, and very 
high ratios, > 1000, of the X-ray to optical fluxes. Their X-ray radiation has been interpreted as 
thermal radiation from cooling NSs. The reason for the lack of observable radio emission remains 
unclear. The simplest explanation is that these NSs are ordinary radio pulsars whose rotational 
and magnetic axes are unfavorably oriented so that the pulsar beam cannot be seen at Earth. 
An alternative explanation is that these objects are not active radio pulsars due to, e.g., too low 
magnetic fields or too slow rotation. Another option is that these NSs belong to the class of 
anomalous X-ray pulsars, also radio-quiet, slowly rotating NSs which apparently have superstrong, 
~ 10 14 ~10 15 G, magnetic fields (see, e.g., Vasisht & Gotthelf 1997). To choose among these 
hypotheses, one should know the rotation period and its derivative, which would allow one to 
estimate the object's age and magnetic field. However, many attempts to find periods in the X-ray 
radiation of these sources have been unsuccessful so far, which indicates that the pulsed fraction 
may be too low to distinguish between a signature of the true periodicity and a noise fluctuation 
in a relatively short observation with a small number of photons collected. On the other hand, 
increasing the observational time strongly increases the number of trial periods (and period 
derivatives in the case of very long observations). This not only makes the search prohibitively 
expensive, but also increases the probability of obtaining large noise fluctuations, which hampers 
detecting the true periodicity. However, if two or more timing observations of a source have been 
carried out, well separated in time, with good time resolution and sufficient amounts of photons 
collected, one may use additional constraints to discriminate between true and false periodicity 
signatures in the separate timing data. We applied such an approach to searching for periodicity 
of the X-ray brightest, radio-quiet NS candidate RX J0822-4300 in the supernova remnant (SNR) 
Puppis A, making use of ROSAT PSPC and HRI observations. 

Puppis A is one of three known oxygen-rich SNRs in the Galaxy. Its kinematic age, 
3.7 ± 0.4 kyr, was estimated from the proper motion of fast-moving filaments in the northeast 
quadrant of the SNR (Winkler et al. 1988). Its distance, d = 2.2 ± 0.3 kpc, was estimated from 
VLA observations in the A21 cm hydrogen line in the direction of Puppis A (Reynoso et al. 1995). 
The central compact source was seen in the Einstein HRI images, located ~ 6' from the kinematical 
center of SNR expansion (Petre et al. 1982). Petre, Becker & Winkler (1996; hereafter PBW96) 
analyzed ROSAT HRI and PSPC observations (exposure times 4 ks and 6 ks, respectively) and 
suggested that this source is an isolated NS born in the supernova explosion. It has no optical 
counterpart to limiting magnitudes B > 25.0, R > 23.6, which implies an X-ray/optical flux ratio 
fx/fB > 5000. No obvious radio counterpart was detected to a 3a upper limit of 0.75 mJy at 1.4 
GHz. Fitting the PSPC spectrum of the source with the blackbody model gives a temperature 
of 0.28 ± 0.10 keV, somewhat higher than 0.10 — 0.18 keV expected from standard NS cooling 
models. The corresponding blackbody radius of ~ 2 km is smaller than the expected radius for a 
NS. No evidence for pulsations with a pulsed fraction larger than 20% was found. 



New observations of Puppis A and its central source with the ROSAT and ASCA have 
been performed recently. We used these data for further investigations of the spectrum and the 
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light curve of the NS candidate. Results of the spectral investigation will be published elsewhere 
(Zavlin, Pavlov, & Triimper 1998b). Briefly, they show that the temperature and radius of the 
NS are compatible with standard models if to assume that the NS surface is covered with a 
hydrogen/helium atmosphere (cf. Zavlin, Pavlov, &i Triimper 1998a). Fitting the observed spectra 
with the NS atmosphere models indicates that the surface magnetic field may be rather high, 
> 4 x 10 12 G, so that one can expect that the NS magnetic poles are hotter than the equator 
(Greenstein & Hartke 1983), which may lead to observable modulation of the thermal X-ray flux. 

2. Observational Data and Timing Analysis 

Among several observations of Puppis A in the public archive at HEASARC/GSFC, two are 
suitable for searching for periodicity of the central source. The first one was carried out with 
the ROSAT PSPC on 1991 April 16 (6.0 ks useful exposure, 53 ks total time span between the 
beginning and the end of the observation). In that observation the NS was off-set by ~ 12'.5. 
The background-subtracted count rate is 0.237 ± 0.006 s" 1 (0.296 ± 0.007 s^ 1 after applying the 
dead-time and vignetting corrections). For the period search, we used 1368 counts extracted in 
the 0.1-2.4 keV range from a 25" radius aperture. 

Another useful observation was carried out with the ROSAT HRI on 1995 November 2-8 
(30.7 ks useful exposure, 580 ks total time span). We evaluated the background-subtracted source 
count rate, 0.070 ± 0.002 s^ 1 (0.080 ± 0.003 after the dead-time correction). An analysis of the 
radial distribution of the source counts does not show significant deviations from the HRI Point 
Spread Function for a point source. To reduce the background contamination, we used 1968 
counts collected from a 10" radius aperture centered at the point source (a 20 oo = 8 h 21 m 57!5, 
5 20 oo = -43°00 / 14"5). 

The ROSAT HRI observations of 1991 and 1992 were too short (2.4 and 1.5 ks, respectively) 
to collect enough point source photons. RX J0822-4300 was also observed with the ASCA GIS 
(1993 July), but a too low time resolution, 65.2 ms. 

To search for periodicity of RX J0822-4300, we converted the spacecraft clock (SCC) arrival 
times of the HRI and PSPC photons to UTC times, and then to barycentric dynamical times 
(TDB). The SCC-UTC conversion was performed with a fourth order polynominal fit of the SCC 
calibration points to UTC. The fit residuals were within a 3 ms range; this corresponds to a scatter 
in pulse phase < 0.1 for frequencies < 30 Hz. 

To perform a global periodicity search of sparse data, tests based on Z^ n statistics (m is the 
number of harmonics included) look most appropriate (Buccheri et al. 1983; De Jager 1994). Since 
the point source spectrum is likely of a thermal origin (PBW96; Zavlin et al. 1998b), the pulse 
profile is expected to be smooth, with main contributions coming from the first harmonic, as 
observed in thermal X-ray radiation from radio pulsars (e.g., Anderson et al. 1993) and predicted 
by NS atmosphere models (Shibanov et al. 1995). This means that the Z\ (Rayleigh) test should 
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be optimal for the periodicity search: the ephemeris parameters / and / are estimated as the 
values that give the largest value of Z 2 = 2N~ l [(^2f =1 cos27r</>i) 2 + (X)£Li sin 2irfa) 2 ], where N is 
the total number of events analyzed, fa = f Aij + / (Aij) 2 /2 is the phase, / and / are the trial 
ephemeris parameters, Aij (i = 1, ...,N) is the event arrival time counted from an epoch of zero 
phase. 

If the arrival times are uniformly distributed (periodic component is absent), then the 
probability density function of Z\ is equal to that of a x 2 with two degrees of freedom, with 
both the mean and the standard deviation equal to 2. The expected (mean) number of peaks 
with Z 2 > Z 2 in J\f statistically independent trials is M exp(— Z 2 /2). This means that if a 
maximum value of Z 2 in an /, / domain is such that a = AAexp(— Z 2 max /2) < 1, periodicity at 
the corresponding / and / can be established at a confidence level of C = (1 — a) x 100%. 

For a sinusoidal signal with frequency /o, frequency derivative /o, and pulsed fraction p, the 
mean and the standard deviation of the peak value Z 2 (/o,/o) can be estimated as Np 2 /2 and 
(2N) 1 ^ 2 p, respectively, provided Np 2 3> 4 and p 2 < 1. This means that the periodicity can be 
established at a high confidence level C only if N and p are so large, and/or M so small, that 
Afexp(—Np 2 /4) < 1 — (C/100%). If the quantity in the left-hand side of this inequality is not 
small enough, the peak corresponding to the true ephemeris parameters may be comparable to, or 
even lower than statistical fluctuations of Z 2 at other values of / and /. In this case, the true 
parameters /o and /o cannot be determined reliably from a single observation. 

The number J\f of statistically independent trials is determined by the domain of /, / values 
chosen for the search, and by the time span T span of a given observation. If | /max - /mini ^ S pan/^ ^ 1> 
then Z 2 is almost independent of /, and one can use a one-dimensional (frequency) grid at a 
fixed / from the range chosen. The number of independent trials is N = (/max — /min)^s pan , and 
the number of grid points should be greater by a factor of ~ 10, in order not to miss the peak. 
In the opposite case, |/ max — / m in|r 2 pan /2 > 1, one has to use a two-dimensional /, / grid, with 
N = (/ max - /min) |/max - /min|r s 3 pan /2 and the number of the grid points of ~ 10 2 AA. 

A reasonable /, / domain for a young, isolated NS is 0.01 < / < 100 Hz, 
— 1 x 10~ 10 < / < Hz s _1 . The 3 ms accuracy of the SCC-UTC conversion hampers 
detection of high-frequency modulation, so that we choose a conservative upper limit, / < 30 Hz. 
For the /, / domain chosen, Z 2 is practically independent of / for the PSPC observation 
(|/|max7 1 span /2 = 0.14), contrary to the HRI observation (|/| max T 2 pan /2 = 17), and the numbers 
of independent trials are M = 1.6 x 10 6 and 2.9 x 10 8 , respectively. Because of the enormous 
number of the two-dimensional grid points needed for the blind periodicity search in the HRI 
data, we have to start from investigating Z 2 for the PSPC events, the first step of our timing 
analysis. We put / = and calculated Z 2 (f) at equally spaced 2 x 10 7 frequencies in the range 
0.01-30 Hz. The result is not encouraging — the highest peak, Z 2 = 32.0, is not much higher than 
several other peaks (Z\ = 30.5, 29.8, 29.1, 27.8), and the probability that Z\ > 32 for uniformly 
distributed arrival times is as large as 0.18. Although any one of these peaks could be caused by 
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real pulsations, with an expected pulsed fraction p = 20.9(-Zf/30) 1 / 2 % for a sine-like pulse profile, 
the corresponding confidence level would be too low (e.g., 82%, or f .3a, for Zf = 32) to claim the 
periodicity is detected. Thus, we have to conclude that only an upper limit on the pulsed fraction, 
~ 20%, can be established from the PSPC data alone, in accordance with PBW96. 

However, with one more set of timing data available, we can examine whether the separate 
PSPC Z\ peaks correspond to HRI Z\ peaks at reasonably shifted frequencies. If a PSPC peak 
at / = /pspc is caused by a periodic signal, then Z\ calculated with the HRI arrival times should 
have a peak at a frequency / H ri = /pspc + fT, where T = 1661 d 00003469 is the time in the TDB 
scale between the epochs (start times) of the PSPC and HRI observations. If such a peak is not 
found for reasonable values of /, we have to conclude that the PSPC peak under examination 
is most likely a statistical fluctuation. Thus, the second step of our analysis is a conditional 
search for periodicity in narrow strips of the /hri 7/ plane. The width of the strip along the /hri 
axis, A/hri = 2/T S p|„ c ~ 4 x 10~ 5 Hz, is determined by the uncertainty of /pspc due to the 
finiteness of the time span. The size of the strip along the / axis is merely |/| max = 1 x 1CT 10 Hz 
s _1 . The number of statistically independent trials for each of the strips can be estimated as 
AT = A/hri |/| max (T s HM) 3 /2^ 370. 

We examined the ten highest PSPC peaks and found no corresponding HRI peaks of a 
high significance for nine of them — heights of the HRI peaks do not exceed Zf = 16.5, which 
corresponds to confidence levels C < 90.4% (< 1.7a). Only for one PSPC peak, Z\ = 30.5 
at /pspc = 13.2875651(8) Hz, a highly significant (C > 99.95%, or > 3.5a) HRI counterpart 
was found, Z\ = 27.3 at / H ri = 13.28378836(7) Hz and / = -2.63(2) x KT 11 Hz s -1 . The 
digits in parentheses indicate uncertainties, q/T spim and 2q/T^ pan , of the last digits of / and 
/, where q ~ (3 ms) /o ~ 0.04 is a phase resolution (cf. Mattox et al. 1996). The values of 
the PSPC and HRI ephemeris parameters are for the (TDB) epochs MJD 48361^50067344 and 
MJD 50022^50070813, respectively. The period and its derivative, for the latter epoch, are 
P = 75.2797300(4) ms, P = 1.49(1) x 10~ 13 s s^ 1 . To verify that the PSPC and HRI Z\ peaks 
are indeed associated with RX J0822-4300 , we applied the Z\ test to comparable numbers of 
off-source events in the PSPC and HRI observations and obtained Z\ values close to the mean 
value expected for a nonpulsed signal. Figure 1 shows the frequency dependences Zf(f) at fixed 
f = fo for the PSPC and HRI data. The main (highest) peaks at / = /o are accompanied by a 
few side peaks at / = fo ± 1.74 x 10~ 4 /c Hz (k = 1, 2, . . .), caused by exposure gaps associated with 
the ROSAT revolution around Earth with the period of 96.0 min. Both the main and the side 
peaks have an additional fine structure with characteristic "periods" l/T S p|^ c = 1.9 x 10 -5 Hz 
and l/T S p^ = 1.7 x 10 -6 Hz. Similar coherent structures can be also seen in the dependence of 
the HRI Z\ on / at a fixed /, with a shortest "period" 2/(T™) 2 = 5.9 x 10~ 12 Hz s _1 . 

Figure 2 shows the pulse profiles extracted by folding the arrival times with the ephemeris 
parameters obtained for the PSPC and HRI data sets. Each of the curves is drawn for the 
corresponding epoch, which causes the apparent phase shift. This shift could be eliminated by 
combining the PSPC and HRI events in one set and determining the ephemeris parameters with a 



-6- 



30 - / =13. 2875651 Hz 

25 1 




•1x10" 



-5x10 



25 
20 
15 
10 

5 



-1x10" 



/ =13. 28378834 Hz 
/ = -2.63xl0" n Hz s" 1 




-5xl0" 4 

f-fo (Hz) 



5xl0" 4 



1x10" 



Fig. 1. — Power spectra around the pulsar frequency /o for the PSPC (1991 April) and HRI (1995 
November) observations, at the frequency derivative / = —2.63 x 1CT 11 Hz s _1 . 
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Fig. 2.- 



Folded light curves for the PSPC and HRI data. 
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much higher accuracy, but the total time span, T span ~ 1.441 x 10 8 s, is so long that this cannot 
be done unequivocally. We verified that there is no statistically significant difference between the 
shapes of the PSPC and HRI pulses. The light curves binned in 12 phase bins have broad maxima 
consistent with a thermal origin of the X-ray flux. The powers of higher harmonics are small in 
comparison with Z\ (3.2 vs. 30.5 and 1.0 vs. 27.3), which demonstrates the pulse profile shape is 
close to a sinusoid and justifies a posteriori our choice of Z\ as a test statistic. To estimate the 
pulsed fractions, we applied the bootstrap method proposed by Swanepoel, de Beer & Loots (1996) 
and obtained p = 22.5 ± 4.2% and 20.1 ± 3.5% for the PSPC and HRI light curves, respectively. 
We also calculated the PSPC pulse profiles at lower and higher photon energies (0.1—1.1 and 
1.1-2.4 keV) and found no statistically significant phase shift or difference of the pulse shapes. 

3. Discussion 

With the period and its derivative evaluated, it is straightforward to estimate the pulsar's 
age, magnetic field, and rotational energy loss. Assuming a pulsar spin-down equation / = —Kf n , 
where n is the braking index and K = const, the characteristic (dynamical) pulsar's age equals 
r = —f/[{n — 1)/] = 16.0/ (ra — 1) kyr, provided that the rotational rate at the NS birth was much 
greater than its current value, /j n 3> /. For n = 3, corresponding to a simple magnetic-dipole 
braking, we obtain r = 8.0 kyr. If we accept that the true pulsar age coincides with the the 
kinematic age of the SNR, 3.7 ± 0.4 kyr, then the braking index should be n = 5.3 ± 0.5. 
To match the pulsar's age and the kinematic SNR age at 2.2 < n < 2.9, observed for the 
youngest radio pulsars PSR B0531+21, B0540-69, B1509-58, f in should be within a 16.8-18.9 Hz 
range (Pj n = 53-59 ms). A remarkably similar situation with the SNR/pulsar age and the 
braking index is found for PSR J1617-5055 near the SNR RCW 103 (Torii et al. 1998; Kaspi et 
al. 1998). The conventional "magnetic field" of the pulsar, derived from the standard formula 
B = 1.0 x 10 12 (PP-is) 1 / 2 G (e. g., Manchester & Taylor 1977), equals 3.4 x 10 12 G for a 
standard moment of inertia , / = 10 45 g cm 2 , and a NS radius of 10 km. The rotational 
energy loss can be estimated as E = —47r 2 Iff = 1.4 x 10 37 erg s -1 . The parameters inferred 
for RX J0822-4300 are close to those of the Vela pulsar, P = 89.3 ms, P = 1.25 x 10" 13 , 
P/{2P) = 11 kyr, B = 3A x 10 12 G, E = 0.7 x 10 37 erg s _1 . 

These results indicate that RX J0822-4300 may be a regular radio pulsar whose radio-quiet 
nature can be explained by an unfavorable orientation of the pulsar beam. Contrary to the radio 
flux, the X-rays observed from RX J0822-4300 originate from the whole NS surface (Zavlin et 
al. 1998b), their pulsation with p ~ 20% may be due to anisotropy of the surface temperature 
distribution and magnetized atmosphere radiation. However, it is expected that the lower is the 
frequency, the broader is the radio beam. Therefore, deep radio observations at low frequencies 
could detect RX J0822-4300, as it has happened recently with Geminga (Kuzmin & Losovsky 
1997; Malofeev & Malov 1997). 

Such a young, energetic pulsar is expected to power a pulsar wind nebula, either in the form 



-9- 



of a plerion or a bow shock. The ROSAT images show two "blobs" , « 2' north and « 1' south 
of the NS, with luminosities L x ~ (1 - 2) x 10 34 (i| erg s 1 in the ROSAT range; the blobs are 
apparently connected to the NS with strips of lower brightness. In addition, the ASCA GIS image 
at higher energies, E > 3 keV, reveals a structure resembling a bow-shock nebula, about 5' — 6' 
east of RX J0822-4300. Although these structures might be manifestations of the pulsar activity, 
one still cannot firmly exclude that they belong to the SNR and are not physically associated 
with the pulsar. Future X-ray and radio observations of RX J0822-4300 and its environment are 
needed to elucidate the nature of these structures. 

It is our pleasure to thank Rob Petre, the referee, for very valuable comments. The data 
preparation for the timing analysis was done with the EXSAS software developed at MPE. This 
work was partially supported through NASA grants NAG5-6907 and NAG5-7017. 
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